THE BELL SYSTEM 

TECHNICAL JOURNAL 

DEVOTED TO THE SCIENTIFIC AND ENGINEERING 
ASPECTS OF ELECTRICAL COMMUNICATION 

Volume 50 November 1971 Number 9 

Copyright © 1971, American Telephone and Telegraph Company. Printed in U.S.A. 

Power Spectra of Multilevel Digital 
Phase -Modulated Signals 

By B. GLANCE 

(Manuscript received May 14, 1971) 

A formula is derived for com pitting the power spectrum of multilevel 
digital phase-modulated signals. The results apply to arbitrary pulse 
shapes and probability distributions providing that the pulses do not 
overlap and are independent. This formula can be applied easily to com- 
pute power spectra of signals phase-modulated by various pulse shapes. 
Several examples are given for rectangular and raised-cosine pulses. 

I. INTRODUCTION 

Power spectra of digital angle-modulated signals can be calculated 
in many ways. The direct way of defining the power spectrum is to 
find the Fourier transform of a sample of the signal on a finite time 
interval T„ . The magnitude square of this Fourier transform is then 
divided by T B and averaged over all possible values of the signal. The 
power spectrum is finally obtained by taking the limit of the previous 
result as T a tends to infinity. Power spectra of binary frequency shift- 
keyed signals have been calculated by this method by W. R. Bennett 
and S. O. Rice. 1 R. R. Anderson and J. Saltz 2 have extended the analysis 
to multilevel digital frequency-modulated signals by using the same 
technique. 

Power spectra of digital phase-modulated signals can also be obtained 
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from the Fourier transform of the signal autocorrelation function. 
This second method has been used by L. Lundquist 3 to calculate 
the power spectra of signals phase-modulated by pulse stream 
]C* a>kg(t — kT). The results obtained in his analysis apply to the case 
of overlapping pulses providing that the random discrete variables 
a k are independent and have identical probability distributions. 

In this paper, a general expression is derived for the power spectrum 
of multilevel digital phase-modulated signals by using the Fourier 
transform technique. The only restriction in these calculations is that 
the signal is modulated by independent non-overlapping pulses. Other- 
wise each level can be characterized by a different arbitrary pulse 
shape and have a different probability distribution. In order to simplify 
the spectral analysis, we show that a multilevel digital phase-modulated 
signal given by 

v(t) = cos M +ZE alg r (t - kT)] (1) 

k t 

can also be written in the case of non-overlapping pulses as 

v(t) =11 aly(t - kT) cos M + PM]. (2) 

k t 

The g T {t) in (1) define the pulse shapes of the different levels. The 
discrete variables a[ can take the values or 1 and are mutually ex- 
clusive in the same time-slot T. In the equivalent expression (2), y(t) 
is a unity rectangular pulse time-limited to one signaling interval T. 
The P r (t) are periodic functions of period T and are respectively equal 
to the pulse shape functions g r (t) in the interval T. 

The spectral formula which is obtained has the same form as the 
expression found by H. C. Van Den Elzen 4 for data signals which can 
be written as a random process: 

<t) = E EdUt- kT). (3) 

k r 

Digital phase-modulated signals written as (2) can be seen as part 
of the class of signals given by (3). The same spectral distribution 
may therefore be expected, once the (jk, r {t) are made explicit. 

The spectral analysis made in this paper follows the method de- 
veloped by Bennett and Rice. 1 A simple compact formula is obtained 
which depends on the pulse shape functions g r (t) and their probability 
distributions. The power spectra are found to be generally the sum of 
a continuous spectrum and of a spectrum made of discrete lines. These 
discrete lines occur at the carrier frequency w c and at frequencies 
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shifted from w c by multiples of the timing frequency. Several examples 
of application are given for 

(i) rectangular pulses of duration I ^ T 
(it) rectangular pulses with finite rise time and decay time 
(Hi) raised-cosine pulses. 

The asymptotic limit of the spectrum as / — > oo is also given for these 
pulse shapes. 

II. LINEAR FORMULATION OF DIGITAL PHASE-MODULATED SIGNALS 

Digital phase-modulation presents some similarity with amplitude 
modulation. For example, binary phase shift-keying by x(t)(ir/2) 
radians is identical to amplitude modulation by x(t) if x(t) is a rec- 
tangular wave which takes the values 1 or —1, thus 

v(t) = cos (u r t + x(t) I) 

= x(t)smw c t. (4) 

The same signal can also be written with T(t) = [1 + x(t)/2] as 

v(t) = r(0 cos (o>,.t + |J + [i - r(0] ^os (« e < - |j- (5) 

In a previous paper, 5 it has been shown that this expression can be 
generalized for any pulse shape which is time-limited to one timing 
period T. Equation (f>) yields in the case of a binary encoding by two 
pulse shapes </,(/) and (/ 2 (0- 

v(t) = r(/) cos [m.i + P.(0] + [i - r(0] cos [ Ue t + p a (t)l (6) 

Pj(t) and P 2 (t) are two periodic functions equal respectively to g^t) 
and g 2 (t) in a timing period T. 

In order to extend this formulation to multilevel digital phase- 
modulated signals, let us consider the multilevel baseband-encoded 
signal 

0(0 - E Zalg r (t-kT). (7) 

A = -oo r = l 

The pulse shape functions g r (l) associated to the different levels can 
be arbitrary providing that 

g r (t) = for > t > T r = 1, 2, • • • , n. (8) 
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The a{ are discrete variables which take the values 1 or at each time 
slot T. The a r k relative to any time slot T are mutually exclusive since 
only one pulse can exist for each time interval T. These conditions 
lead to 

£ al - 1 for kT ^ t£ (fc + 1)T 

T- 1 

Jfc = _ „ f . . . 0, • • • + oo . (9) 

For a signal made of independent pulses, the a r k relative to different 
time slots form a discrete stationary random process. 

The baseband signal (7) can be seen as n mutually exclusive signals, 
each one related to a particular level as shown in Fig. 1. For example 



9>(t)=zza k g r (t-kT) 
k r 
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Fig. 1 — Multilevel baseband-encoded signal. 
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the level r is given by 



#r(0 - E «^r(i " WO- 



(10) 



Each <f> T (t) can itself be seen in Fig. 2 as the product of a unity rec- 
tangular wave by a periodic function, thus 



*,(*) =PXt) !Laly(t-kT). 



(11) 



P r (t) is a periodic function of period T equal to g T (t) in the interval T. 
y(t) is a unity rectangular pulse time-limited to one interval T. Let 
us set 



r r (0 = Zaly(t-kT), 

k 

the baseband signal becomes 

4>(t) = Er r (fl-p r (0. 



(12) 



(13) 



A signal phase-modulated by an n-level digital signal such as (7) 
can be written from (13) as follows: 



v(t) = coskt+ Zr r (/)-P,(01. 



(14) 



The n unity rectangular signals T T (t) as given by (9) are mutually 
exclusive for all values of t: 



Zr r (0 = i. 



(15) 



<p r (t)= s a£g r (t-kT) 

K 



/S, ,/s,/\ 



£v 



PpCt) 



r P (t)=sa^ y (t-kT) 



Ai/N 



/Si /Si 



T 2T 3T 4T 5T 6T 7T 8T 9T 10T t 

Fig. 2 — Baseband signal corresponding to the level r. 



2862 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1971 

As a result, (14) is also equal to 

v(t) = E r r («) cos M + P r (t)l (16) 

The modulated signal written in this form presents some similarity 
with amplitude modulation. Note that the argument is independent 
of the random variables a r k . This equivalent formulation simplifies 
greatly the spectral analysis of multilevel phase-modulated signals. 

III. POWER SPECTRA OF MULTILEVEL PHASE-MODULATED SIGNALS 

3.1 General Considerations 

The spectrum calculations are made under the following conditions: 

(i) The baseband-encoded signal is made of independent pulses. 
(ii) The levels are defined by an ensemble of functions g r (t) which 
are time-limited to an interval shorter than or equal to the 
timing period T. 

Otherwise the g r (t) are arbitrary and can have different probability 
distributions. Let p r be the probabilitj' of occurrence of the level r. 
From (10) and (12) 

Pr = <«[.), (17) 

and from (15) 

EPr-1. (18) 

The power spectra, calculated by the Fourier transform method, 
are given by 

<7(/) = Lim ^= <| X(/, NT) | 2 >. (19) 

NT— oo *» L 

X(f, NT) is the Fourier transform of a sample of the signal in the 
interval NT. (\ X(f, NT) | 2 ) is the expectation taken over all possible 
values of the signal in the interval NT. 

The Fourier transform X(f, NT) is calculated by using the linear 
formulation (16), thus 

X(f, NT) = E f r r (0 cos M + P r (t)]e~ ia ' dt. (20) 

r •'0 

X(f, NT) can be separated in two parts which correspond respectively 
to the positive and negative frequencies: 
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X(f + ,NT)-lE (^ ^ r (0• ,IC -- ,,+ ' rC,), dt (21) 

r Jo 

and 

X(/_ , OT) = | £ r' r r (^-'' ,(Uc+u)t+Pr<m eB. (22) 

r •'0 

Note that X(/_ , W) = X * (-/+ , OT) [* indicates that ; is made 
-j in (22)]. 

The Fourier transform (22) relative to the positive frequencies 
can be written from (12) as 

X(J + ,NT) = \ E t al r '«"<--»♦"'«>' dt. (23) 

* = r=l JkT 

Let us set the new variable of integration y = t — kT. From the perio- 
dicity of the functions P r (t) 

P r (y + kT) = P{y), (24) 

which gives 

X(j + ,NT) = \ £ Z atf"'""" ' f • "•—>•♦*<"' ty. (2.5) 

«.- = () r=l •'0 

In the limit of integration P r (y) = g r (y); g r (y) can then be substituted 
to P r (y), thus 

X(j + , NT) -i'fi; atf*""-" /" e '«— «>»♦"<»» ,/,„. (26) 
Let us set 

F r (/ + ) =^ fe'-"-"^""^ (27) 

and 

F r (f_) = \ f'e-'ic.+.^wi^ (28) 

.s Jo 

The Fourier transform (20) takes the form 

X(j, NT) = E I) aSfe"' 1 —" 'F r (/ + ) + c - / " ,( - + - , F r (/_)} . (29) 

k-0 r=l 

3.2 Power Spectrum Calculations 

The power spectrum calculated from (19) is the sum of four terms 
given by 
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j /A-i n 

?(/) = Lim rj— { 2 I] alal. 



oty 



F r (/ + )-F*(/ + y m -' ,(w '- u) 
4-^ r (/-)-^(/-)c"'' T(fc "' ,(lJf+w) 

+F r (/ + )-Ff(/_)e'" :7 ' (0 ' c "" ) + ' 7 ' ( ' Je + ' J) ' 
[4-^ , r (/_)-F,*(/ + )e _yu ' 7 ' (Uf + u) + ' r(Uf " <J " 



(30) 

The third and fourth terms give the foldover between the positive 
and negative frequency parts of the spectrum. This effect is negligible 
when the bandwidth spectrum is much smaller than the carrier frequency 
which is the case in most applications. The spectral distributions are 
then accurately described by the first and the second terms of (30) 
which give respectively the power spectra relative to the positive 
and negative frequencies. The power spectrum relative to the positive 
frequencies yields in that case 

f/(/ + ) = Lim ^ (Z p aI.a; e ' ru '-' ,(w '- u, F r (/ + )^(/ + ))- (31) 



The ensemble average is taken on the random variables a r k a", . The 
bracket sign can therefore be introduced on the summation as 

GQ + ) = Lim Jl. jj) 2 WWy* c *- ,H —" X(/ + )-TO + )}- (32) 



(33) 



From the definition given in Section II: 

\p r p, if k ?* I and r ^ s 

. . ,. pi if k 9^ I and r = s 
(a k a,) = -T 

p r if k = I and r — s 

.0 if k = I and r ^ s 

where p T and p, are the probability distributions of the levels r and s. 
After substitution of the terms (ala',) by their values given by (33), 
a summation is made on the variable I. The result gives 



G'(/ + ) = Lim I £ |y Vp^ S' r + y p r p. 



•z 



(AT - fc) 

N 



cos kT(a c - o>))F r (f+)F*(j + ) 



jFr(f + )F*(f*)\ , 



(34) 
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5; is the Kronecker symbol given by 

1 if r = s 



K = 



.0 if r ?± s 



Since the limit operation acts only on the second term of (34), G(/ + ) 
can be rewritten as 

G(j+) = |_E {y Vp,p. K + 2p rPl Lim ^ 

= 1 

. g (JV - fc) cos kT(a e - co)JF r (/ + )F.(/ + )]. 

Summation on the discrete variable fc yields for the series on A; 



(35) 



Lim t|p E (N - k) cos />T( Wr - to) 



-V T— » 



ot 



= Lim t= 



-N + 



. , .ATT 

5 in (cu r — oj) -77- 



sin" (co c — «) — 



(36) 



s «(/. - / + f 



forms a set of Dirac functions which give a series of spikes: 



E Af. - / + 



m 



m 



+ « if /-/. + f 



(37) 



m 



if /f/. + f 



The substitution of (36) into (35) yields finally for the power spectrum 
relative to the positive frequencies 



G(j + ) =f\EPr\FrV + ) I 2 - 



E PrF r (/ + ) 



1 +0 ° 
+ ?-£ 



E PrF r (/- 



K/.-/ + ? 



(38) 



The first term of (38) gives the continuous power spectrum. The second 
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part gives the spectral discrete lines. These lines occur at the frequencies 



m 



/ = U + f (m = 0, ±1, ±2, 



) 



if 



Z VrFrii.) 



m 



^ at / = U + f 



The final formula (37) depends on two sets of parameters: 

(i) the probability distributions p r of each level 
(it) the pulse shape functions g r (t) of each level from the Fourier 
transforms (27). 

When the n levels have identical probability distributions, (38) 
becomes 



(m = \^±\F r (U)?-^ 



^ T 2 J? m n 2 



E Fr(f+) 



E Fr(f + ) 



i , , m 



(39) 



The pulse shape, in most cases, is given by a unique function g(t) 
for all the levels. Each level is then denned by a different pulse amplitude. 
In this case, g(t) can only have one maximum value which is assumed 
to be equal to one. The levels are then given by the ensemble a r g(t), 
r = 1, 2, • • • ,11, where the a T are equal to the peak phase deviations. 
Several examples of power spectra will be calculated for this important 
case. 

IV. EXAMPLE OF POWER SPECTRA 

4.1 Rectangular Pulse Shape of Duration T 

The application of (38) to a signal phase-modulated by rectangular 
pulses of duration equal to the timing period T is straightforward. 
The Fourier transforms (27) are given in this case by 



F r (U) = \ f e m , —' 1 dy 

- Jo 



J(u e -u)(T/2) 



sin (w c — ui) — 



(40) 



where the a T are the phase deviations of the different levels. Equation 
(40) applied to (38) gives 
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T 



sin (w f — <ji) — 
(oj r — w) — 



1 - 



£p,e"" 



Ep/"' 



*(/. - /)■ (41) 



The spectrum has only one line at the carrier frequency. 

If one assumes that the levels are equidistant and equiprobable 
such that 



, r- 1 
a T = a„ -f- a 



Vr = 



one obtains 



Spy 



and the spectrum is given by 

Tl 2 



r = 1 2 



n . a 

sin — 
2n 



n 



(42) 



(43) 



0<M = 



7' 



sin (co c - w) - 



(«. 



")o 



1 - 


n 2 


sm - 


■ + * 


sin - 


a 

sin — 
2/i 


. a 

sin — 
2n 



S(f e -f). (44) 



In polar modulation, the a r can take positive and negative values. 
The term S"-i P' e '°' can tnen become equal to zero for an infinite 
number of solutions. All these solutions give the same power spectrum, 
but without the discrete line: 



T 



i 

sin (co r — oj) — 

(co r — w) 5" 



(45) 



For example, the polar systems made of q pairs of equidistant and 
equiprobable levels such that 
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^2 e' a ' = ^2 cos a T = 
are given by 



2 levels polar modulation <a ± , = ±- 



4 levels polar modulation 



6 levels polar modulation 



a ±1 = 


-I 


a ±2 = 


4 


a±i = 


±5 



a ±3 = ± J 



etc. 



All these modes have the same spectrum given by (45). 

4.2 Rectangular Pulse of Duration t < T 

Let us now consider a rectangular pulse shape of duration t < T. 
Assume the pulses are located symmetrically in the interval T as 
shown in Fig. 3. The Fourier transforms F r (f+) are equal in this case to 

, f r (T-r)/2 



+ •*-' / e i{ "<-^ dy + e iu "-" u dy\- 



(46) 



g(t) 




Fig. 3— Rectangular pulse with t < T. 
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The result of (46) applied to (38) yields for the continuous part of 
the spectrum 

T 



Ge(j + ) =l\HVr 



sin (w r — w) — sin (w r — w) r 



r 

(«« — w) « 

T 



(«„ - w) - 



(1 - «"") 



sin (co r -co) — sin (co r — w) r 



(w« — «) o 



The discrete lines are given by 



(w e - w) - 



i-Ep/ 



(47) 



= 1 g 



sin mir 
mir 



Sin W?7T TTi 
T_ i_ 

T " r 

mir j, 



1- EPrC 



•1/.-/ + 



(48) 



There are always discrete lines in this case. As an example let us con- 
sider the case of a polar phase modulation with two equiprobable levels 
given by a r = ±ir/2 radians. Equations (47) and (48) give in this case: 



w.)-ffe 



sin (w r — a>) - 

(»c — ") - 



| +00 

+ i 2 

* m = — r*i 



si n mir 
mir 



Sill W7T ^ 

T_ i_ 

mir ji 



l/ e -/ + | 



(49) 



For 7 1 — t « T the continuous spectrum is practically the same as 
the one obtained for pulses of duration T, except for the discrete lines 
at / = f r + m/T (m = ±1, ±2, • • ■ ). In both cases the spectrum 
decreases as l//~ as f —* x . 

4.3 Rectangular Pulse Shape with Finite Rise Time and Decay Time 
We consider in this example a rectangular pulse shape of duration T 
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get) '■ T = T " 2S 



s 



T-ST { 



Fig. 4 — Rectangular pulse with finite rise time and decay time. 



with a finite rise time and decay time as shown in Fig. 4. The calcula- 
tions are made for the case of two equiprobable levels with peak phase 
deviations equal to ±r/2 radians. We assume that the rise time and 
decay time are both equal to s = (T — r)/2; r is the top pulse length 
The Fourier transforms F r (f + ) are given in this case by 

• T-» 



Fl(/+) = 2 {/„' a" 1 -- **"*"" dy + * n * m \[ 



dy 



(50) 



F 2 (f+) can be obtained from (50) by changing the sign of r/2. The 
calculations yield for the continuous power spectrum: 



<ww-ili 



The discrete lines are given by 
1 ItY: 



sill (<a e - co) - 
L 1 _i_ 

T ( ^ r 


^ COS (<J) r - 


-«) 2 

> 


1 («r — ") g 


2 - (i) 1 





(51) 



a,(/ ( ) = 



4 \4/ VT 



cos" [ 2irm ™ 



[(- 



T 



*/.-/ + ? 



(52) 



The power spectrum decreases in this case as l// 4 as / — » «> . Figure 5 
shows the spectrum given by (51) and (52) which is calculated for 
s/T = 0.1. 

4.4 Raised-Cosine Pulse Shape 

In the case of raised-cosine pulse shapes, the spectra are computed 
for polar modulation. We assume that the pulse duration is equal to 
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the timing period T. The pulse, of maximum height 1, is then given 
by g(t) = sin 2 fl^/2 with ft = 2v/T. A first example is calculated for 
a two-level system with different probability distributions. The second 
example is calculated for a multilevel system with equiprobable level 
distributions. The two-level and four-level cases are then compared 
for peak phase distributions which give ^" e'"' = 0. 

4.4.1 Two-Level Polar Phase Modulation with Raised-Cosine Pidses 

Let ±2a radians be the peak phase deviations of the two levels. 
The Fourier transforms F r (/ + ) are given in this case by 



£i Jo 



ilu/-2] 



dij . 



(53) 



F 2 (f + ) can be obtained from (53) by changing the sign of a. 
The term 2ja sin 2 Q.ij/2 of (53) expanded in Bessel series gives 



0.25 




f c+ l f c+ f f c + ^ 



Fig. 6 — Power spectrum for two-level polar phase modulation with rectangular 
pulses with finite rise time and decay time. 
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e 2i *' inta " 2 = e"\j.(a) + 2 2 (-l)V 2n («) cos 2nQy + 2j 



• £ (-l)V 2n _,(a) cos (2n - 1)Q|^- (54) 

Substitution of e 2 '° sin ' fh "' 2 by its Bessel series expansion (54) into (53) 
yields 



F(f + ,±a) = |><— -><™-> 



sin (oj c — co) — 



T 

(«« — w ) ^~ 



±2j f:(-l)V 2n _,(a) 



.- co) 2 - (2nJ2) 
(co c — co) 



(55) 



, (co c - co) 2 - (2n - l) 2 fi 2 

Equation (55) applied to (38) yields for the continuous part of the 
power spectrum: 



O c (U) = PilhT 



T 

sin (co r — co) — 

T~ 
(&>„ - co) - 



Ma) + 2 Z (-l)V 2B (a) 

+ 2 Z(-l)V 2 „_,(a) 



(co c — co) 



(co r -co)' - ( 
(co c — co) 



2^2fJ Sina 



Noting that 
and 



(co e - co) 2 - [(2w - 1)Q] 2 



J„(a) + 2 £(-l)V 2n (a) = cos a 



o cos a 



(50) 



(57) 



2 E(-l)V 2n _,(a) = -sin a, 
(5(5) can be rewritten as 



G e (j + ) = 4 Pl p 2 T 



T 
sin (co r — co) — 

(co r — «) =■ 



E(-DV.(a) 



sin (or + — 



1 " ' C ^f 



(58) 
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The discrete lines are given by 

M+) = \ m £ M J-(«)\} - 4P.P2 sin 2 [a + m |) *(/, - / + f )■ 

(59) 



The asymptotic limit of (?„(/+) as / 

(58) for | u e — co | > ?j?0, thus 



oo is obtained by summing 



Lim G e (j + ) = 4p ]P /r< 



T 

irasin (co r — co) — 



(w fl - co) —J 



(60) 



which shows that G„(J+) ~ l/f as / — » ». 

The power sharing between f/ r (/ + ) and Gj(/ + ) is a function of a and 
]hp 2 • Since pi + p 2 = 1, G>(/+) = if pi or p 2 = 1 as expected. All 
the power is then contained in the discrete lines (case of a signal phase- 
modulated by a periodic function). 

The condition p t = p 2 = \ (identical probability distributions) 
gives a maximum for the continuous distributed power: 



and 



P e (/ + )=| ra L/^)sin 2 (« + m|) 
P*(f + ) =| E •/;,(«) cos 2 (a + m|)- 



(61) 



The power becomes equally divided between both parts of the spectrum 
for a peak phase deviation of ±t/2 radians and pi = p 2 ■ |. For 
this particular case, 8 the series expansion of (55) and (59) can be limited 
to the first two terms. A good approximation for the power spectrum 
is then obtained bv 



T 



sin (co r — co) -^ 



T 

(co, — co) — 



J, \l) ~ 2 ( w ' " l; 



'4i 



>J& 



_(w e — co)" — 12" 



; r - co)" - (20) _ 



and 



Figure 6 shows the power spectrum given by (G2) and (63). 



(62) 



(63) 
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Fig. 6 — Power spectrum for two-level polar phase modulation with raised-cosine 
pulses. 



4.4.2 Multilevel Polar Phase Modulation with Raised-Cosine Pulses 

The results obtained in the previous section can be generalized in 
the case of multilevel polar phase modulation with identical probability 
distributions. Let ±2a r radians be the peak phase deviation of a multi- 
level system of q pairs of levels (r = 1, 2, • • ■ , q). The Fourier trans- 
forms from (55) are equal to 

T 
T sin {oi c — w) - 
F r (f + , ±« f ) = le*""— >«v«*-,i __ _^_ {Ar ±jBrl m 



(w« - w) 2 



with 
and 



A, = J {a r ) + 2 £ (-l)V 2B (« f ) far, m) ' 

i (« c -co) — (Jnii) 



5 r = 2Z(-DV 2 „_,(a r ) 



(to,. — 0))' 



(co c -co)" - l(2n - 1)12]" 



(65) 
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The continuous power spectrum calculated from (38) with (64) and 
(65) yields 



T 



sin (w c -a) g 



•ft £ [AJ + «] - A [Z (A, cosa, - B r sin« r )]y 

into account (57), A T and B 
A r = cosa, + 2 £(-l) n 



Taking into account (57), A T and Z? r can be rewritten as 



1 - 



and 



u r — (J) 



B r = -sin a, + 2 £ (-1)" 



2nfi 
J 2 „M 



(2w - l)flj 



1 - 



Substituting of ,4 r and fi r by (67) into (66) yields 



G.(f + ) = 7 1 



v 



sin (aj r - w) g 
(u e - a) - 

E(-l)V ta (Or) 



1 " =V-^ 



E(-Dv 2 „-,(« r ) 



[ (tf f — a) " i _ I a,f ~ " I 

2nQ J L(2n - 1)0 J 



Z (- l)V 2 „(a r ) cos a r S (- l)V 2n -,(a,) sin a r 



w r — w 



2wQ 



(66) 



(67) 



(68) 



L(2w - 1) 
The discrete lines, obtained from the second term of (68), are given by 

O d {U) = "^ m E { Z JJflr) cos (a, + m |)} 2 fi(/ e - / + f ) • (69) 



4<f 

For a four-level system given by 



2a, = ±- 



2a-. = ± 



3* 
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Fig. 7 — Power .spectrum for four-level polar phase modulation with raised-cosine 
pulses. 



the series expansions in (6S) and (69) can be limited to the second 
term. A good approximation of the power spectrum is then given by 



T 



T 

sin (co r — to) — 

(<a c — w) a 



J' 


©+' 


■(1) 


\b 


\ 12 


<' 



41 j + ./<| 



1 - 



co r — CO 
212 



./ 



/ T 1 "" i r I •^ 7r 1 3t , (w\ . T . r (St\ . 3t 

AS/ 8 2 \ 8/ C0S 8" ' ° ' Sm " ' "°~ ' Sm "° 



1 - 



co, — OJ 

212 



8/ 8 



1 - 



8/ 8 



C0 r — CO 



(70) 
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and 



<!,l(f+) ~ 16 m ? a \ Jm \l) coa (i + m l) 



+ J B ,(|)cos(|+m|)} 2 5 (/„-/ + f). (71) 

The spectrum given by (70) and (71) is shown in Fig. 7. Note that this 
spectrum is almost identical to the spectrum of the two-level system 
calculated with 2a r = ±ir/2 and pi = p 2 , shown in Fig. 6; in both 
cases 2r-i e ±2 ' ar = 0. This result presents an analogy with the cal- 
culations made for rectangular pulses of duration T w'hich satisfy 
the same condition. 

v. conclusion 

A formula has been derived for computing the power spectrum of 
multilevel digital phase-modulated signals. The results apply to arbi- 
trary pulse shapes and probability distributions providing that the 
pulses do not overlap and are independent. 

This formula is easy to apply to various pulse shapes. Several examples 
are given such as: 

(/) rectangular pulses of duration smaller than or equal to the 

timing period 
(it) rectangular pulses with finite rise time and decay time 
(Hi) raised-cosine pulses. 

Approximate results can be obtained for more complex pulse shapes 
by the following method. The spectrum is a function of the pulse shapes 
g T (t) by the Fourier transforms 

F r (j) = I /"V ,lM < +, " (0, e-' 2 "d/- 

T can be segmented into a large number of intervals in which g T (t) 
is approximated to its average value in this interval. F T (j) is then given 
by the following series expansion 

T 

F T {j)c*J- ^-»nrn„ =?.£«»' 

-ii , . l t-n 

(co r — co) — 
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with 

b k = r 
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